DNA denaturation and wetting in the presence of disorder 
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We present a precise equivalence of the Lifson-Poland-Scheraga model with wetting models. Mak- 
ing use of a representation of the former model in terms of random matrices, we obtain, in the limit 
of weak disorder, a mean-field approximation, that shows a change of the critical behavior due to 
disorder. 
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Quenched disorder plays a crucial role in phase transi- 
tions. It can either modify the nature or the order of the 
phase transition [ll-Qor affect the value of the critical 
exponents 0-0] ■ I n this letter we consider the Lifson- 
Poland-Scheraga (LPS) model of DNA denaturation, in 
the limit of weak quenched disorder. In the ordered case, 
the order of the phase transition depends crucially on a 
parameter a @. It is continuous if 1 < a < 2, while 
the transition turns to first order if a > 2. Recently, it 
has been rigorously proven Q that the disordered model 
exhibits a continuos transition for all values of a. The 
values of the critical exponents are unknown, except in 
the case 1 < a < 3/2, where they coincide with the expo- 
nents of the ordered case [l(|, in agreement with Harris' 



criterion . On the other hand, Giacomin and Toninelli 
flpj ] have pointed out the equivalence of the LPS model 
with wetting models. The standard disordered wetting 
model 0, [f| has been thoroughly studied and its solu- 
tion is still challenging. In this letter we establish a 
precise connection between a Levy type wetting model 
and the LPS model, so that a first order transition is 
possible in the ordered case. Our results are based on 
a representation of the disordered LPS model as a ran- 
dom matrix model, making contact with a localization 
problem. This representation allows to obtain an exact 
supersymmetric formulation of the model, although here 
we discuss a less rigorous one, based on replicas. We 
obtain a mean- field approximation, which confirms that 
the critical exponents are unchanged ifl < a < 3/2, 
but if 3/2 < a < 2 the exponents are independent of 
a. However, if 2 < a < 3 the phase transition remains 
first-order, thus indicating that in this case fluctuations 
play a crucial role. 

The disordered case - DNA denaturation is described 
in the LPS model by the presence or the absence of an 
H(Hydrogen)-bond at the lattice site x in a chain of 
length N. A weight X x is associated to the H-bond at 
site x, whereas the absence of H-bonds in a sequence 
of sites of length j is weighted by the algebraic factor 
a(j) = cj~ a . In fact, this model is equivalent to a wet- 
ting model Q. Indeed, let us consider a random walk 
that is constrained to the upper half-plane, with a prob- 
ability p(m) of making a jump of size m. We can assign 



a weight X x at any time x the random walk touches the 
horizontal axis. There is a probability a(J) that during 
time j the walk does not touch the horizontal axis. For 
the usual random walk, it is well known that asymptot- 
ically a(j) = cj~ 3 / 2 . More generally, we can prove that 
the same results holds if J2 m P( m ) m2 < 00 • On the other 
hand, if the walk is of Levy type, i.e. if asymptotically 
p(m) — pm~ a with 1 < a < 3, then one finds asymptoti- 
cally a(j) = cj^^ 1 . This bridges the Levy type random 
walk with the LPS through the relation a = > 3/2. 
In both cases the thermodynamics is determined by a 
partition function q(N), that is given by the sum over 
all possible, properly weighted configurations. A crucial 
quantity in the wetting problem is the density of contacts 
with the horizontal axis, p, which corresponds to the den- 
sity of H-bonds in the LPS model. The phase transition 
is characterized by the fact that the order parameter p 
vanishes beyond a critical value of the control parameter. 

In both models, the partition function q(N) obeys the 
recursion formula 



N-2 



\ N q{N) = J2q(N-l- x)a(x) + 6 



N,l 



(1) 
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where for physical reasons A^ > and a(0) = 1 in 
the LPS model. The ordered case corresponds to X x = 
A , Vcc: it can be easily solved by computing the Laplace 
transform of q(N). One can then deduce from it the 
pressure p = liniAr^oo log q(N) and the density p. One 
finds that p — if A > A c = Y^=i a ( x )- If A < A c , near 
the critical point the density vanishes as p = p[X c —X) a ~ 1 
if 1 < a < 2. On the other hand, for a > 2 one finds 
p{X c ) > 0, so that the transition is first order with a 
jump in the order parameter. It is worth stressing that 
in the Levy type random walk one has a first order tran- 
sition for I < a < 2. In we have also computed the 
structure factor of the LPS model. We found that a basic 
length scale £ is associated to the critical behavior of the 
model, 




if 1 < a < 2 
if 2 < a < 3 
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In the limit N — > oo we find that the structure factor, 
S(k), has a scaling limit 



lim 

fc->0,£-> 



S(k) 
5(0) 
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where -F Q is a universal function that depends only on the 
scaling exponent a, thus showing that a is a universal 
critical exponent associated to this class of models. 

The disordered case - According to some recent contri- 
butions by Giacomin and Toninelli [j| the disordered ver- 
sion of the LPS model is expected to exhibit a continuos 
transition over the whole range of values of a. In order 
to tackle this problem we exploit another useful repre- 
sentation of the partition function as a determinant of a 
random matrix. Indeed, let us define a perturbed parti- 
tion function, q(N,e), where we replace in the recursion 
formula (JTJ) a(x) by (1 + e)a(x). This allows to obtain 
the average density p(N) as the logarithmic derivative of 
q(N, e) with respect to e in the limit of vanishing e. Let 
sus ntroducc the N x N matrix 



C € N {x, y) = (1 + e)a(x - y) - \ y 5 y ^ + i, 



(4) 



where (x, y) € [1, N] and we have assumed that a(x) = 
if x < — 1. It can be shown that the partition function 
can be rewritten as 



N 



q(N,e) = ]JX- 1 detC% 



(5) 



Then the pressure and the density at finite volume can 
be expressed as follows 



P N = ±(\nq(N,0)) 

PN = ^< ln ^e))|e=o 



(6) 
(-) 



where the symbol (•) indicates the average over disorder. 
Let us observe that if \ x > ,V z then det Cjy is also 
positive. 

On the other hand, one can extend these formulae to 
the case where is a gaussian variable with a given 
average (A) = A and with a given variance v = (A 2 .) — A 2 . 
For this purpose, we use the equality 



detC^ = [(detC^'Cjv] 



(8) 



In this way the model can be generalized to gaussian dis- 
order by redefining the finite volume pressure and density 
as follows 



P 



N 
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N 



— F(lnA 2 ) 
2N ^ x xl 



PN = ^^(lng(iV,e))| £=0 



(9) 
(10) 



where 



D N = ^(lndet(C^) t C^) 



(11) 



We want to point out that the quantity H = (C%YC% 
can be considered as a random Hamiltonian. Solving 
the problem of the thermodynamics of the disordered 
LPS model amounts to compute the density of states of 
this random Hamiltonian. In this context, the denatura- 
tion transition we are looking for can be interpreted as 
a passage from a phase of localized eigenstates (p > 0) 
to a phase of extended eigenstates (p = 0). This is a 
priori surprising because in this case the properties of 
the density of states reflects the localization transition, 
amounts to a localization transition from localized to ex- 
tended eigenstates, corresponding to p > to p = 0, re- 
spectively. Usually, in this class of random Hamiltonian 
models the average over disorder can be performed by 
two main techniques: the supersymmctric method and 
the replica trick. The former method is certainly more 
rigorous and we shall present it in a future publication. 
In what follows we shall use the method of replicas, be- 
cause it yields a more transparent interpretation of the 
results. In practice, we aim at computing the quantity 



J N 



((det(ChYCh) n ) 
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so that we can obtain an explicit expression of the nor- 
malized pressure through the relation 



P C N = lim — In 



>o 2n 



N 



(13) 



On the other hand, Z e N can be rewritten as follows 

Z<k = J PVKexp [i>v>4*\ C e N (x, y))> (14) 



where 



2n 



[tpy^x] = X! V'y("T.)'0x(m) 



(15) 



m— 1 



where {ip y (m) and ipxijn)} are grassman variables. Now 
we can consider X x as gaussian random variables or per- 
form the weak-disorder limit. One obtains the expression 

Z% = /X^exp(E Xll , D%{x,y)) 

xexpdEjV^+i] 2 ) (16) 

where 

P>n(x, y) = 0- + e)a(x - y) - A^,x+i (17) 
Making use of the identity 

exp(-|trifj)= (18) 
2 J dQ x exp (-trQ 2 - i^/vtrQ x H x ) , 
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where H x and Q x are symmetric (m x to) matrices and d 
is a suitable normalization constant, we can linearize the 
term trifj = — [V'xjV'x+i] m Eq. (fl"6f and thus 
eliminate any dependence on the grassman variables. We 
finally obtain: 

Z e N = w=r I Ell"! 1 dO* exp(- E x ^Ql 

det [D^ + 2i^/vK]) (19) 



D e N (x,x') and 



Here D e N (x,m;x' ,m') = 
Kjsr(x, m; a;', to') = (Jj./ ia . + i(5 x (m, m'). Since we want to 
focus our attention to the weakly disordered case, we 
can perform the mean-field theory of the problem at 
hand that can be obtained by the following steps: 

• (si) one has to first shift the random matrix: Q x — > 
Qx + i\[vq x \ 

• (s2) the determinant in Eq. (fH?l) can be rewrit- 
ten as dctDjj x exp(tr(ln(I + 2i y /v(D t N )- 1 K)) 
where D e N (x, to; x', to') = D € N (x,x')8 m ,m' - 
2vq x 8 x > , x+ \5 rn.m' 

• (s3) one can finally expand the logarithm to first 
order in v and choose q x in such a way that any 
linear term in Q x disappears from the expression 
of Z% (see Eq.(fT9f). 

One can decompose the partition function into one part 
corresponding to the mean-field contribution (MF) and 
another one corresponding to the fluctuation one (Fl), 
namely 



7 e _ rra MF 7< l Fl 
— N 



where 



MF 



J N 



= det D% exp 



2nv i 



(20) 



(21) 



and 



Zh Fl = l^r J IL dQ x exp(- £~ 2 E {xu ... ,,„} 

s n (xi,--- ,x n )txQ Xl ■■■Q Xn ) (22) 

with 

s 2 {xi,x 2 ) = S XuX2 - 2v{D e N y i (xi 1 x 2 ){D e N y 1 {x 2l xi) 

s n (x!,-- - ,x n ) = i(-2^)™nj=i(^ ) w) _1 (^ J ^+i) 

(23) 

and the latter expression holds for n > 3. Let us il- 
lustrate the mean-field results we obtain after perform- 
ing the limit N — > oo, in such a way that q x — > q and 
A = X + 2vq. The resulting self-consistent equation reads 



d0 
2^ 



cij exp i6(j + 1) — A 

3=0 



(24) 



and the corresponding density is 
p= -1 + Xq 



(25) 



Accordingly, we find that the critical point is shifted by 
disorder 

A c = A c + (l-p c )^ , Co c = 0fora<2) (26) 

and, consistently with the Harris criterion (see [HI), we 
find that the critical exponent is the same as in the or- 
dered case for l<a<3/2, in formulae 



P = p(A c - A) 



On the other hand, 



for | < a < 2 we find 



2v 



A) 



(27) 



(28) 



i.e. the critical exponent is 1 all over this range of values 
of a. On the other hand the Mean-Field solution still 
predicts a first-order transition for 2 < a < 3. This in- 
dicates that higher order correction to the MF solution 
have to be taken into account in order to confirm that the 
transition becomes continuous also in this range in the 
presence of disorder. We guess that an effective technique 
to tackle this problem should be based on a renormaliza- 
tion approach, directly applied to the partition function 
of the disordered model. Since we are dealing with a ID 
models we expect that a decimation procedure can be 
worked out. This notwithstanding, the functional forms 
that appear in the above equations seems to challenge 
the possibility of a fully analytic, even if approximated, 
study. 

Conclusions. Our approach offers the possibility to 
compute exactly the critical exponents if we take into 
account fluctuations in the weak disorder limit. The case 
a > 2 is particularly challenging, since the rigorous work 
of [H shows that large deviations, in the probabilistic 
sense, are responsible for the continuous nature of the 
transition. 

Acknowledgments. RL acknowledges financial support 
from the italian MIUR-PRIN project n. 20083C8XFZ. 



[1] Y. Imry and S.-K. Ma, Phys. Rev. Lett. 35, 1399 (1975). 
[2] M. Aizenman and J. Wehr, Phys. Rev. Lett. 62, 2503 

(1989); Commun. Math. Phys. 130, 489 (1990). 
[3] C. E. Soteros and S. G. Whittington, J. Phys. A 37, R279 

(2004). 

[4] G. Forgacs, J. M. Luck, Th. M. Nieuwenhuizen, and H. 

Orland, Phys. Rev. Lett. 57, 2184 (1986). 
[5] B. Derrida, V. Hakim, and J. Vannimenius, J. Stat. Phys. 

66, 1189 (1992). 
[6] G.Giacomin, "Disorder and critical phenomena through 

basic probability models", in Ecole d'ete' de St-Flour 

XL2010, 2011 Springer Verlag. 



4 



[7] F. Caravenna, F. den Hollander and N. Petrelis, "Lec- 
tures on random polymers", arXiv:1106.6204vl(math 
PR)30 june 2011. 

[8] Y. Kafri, D. Mukamel, and L. Peliti, Phys. Rev. Lett. 85, 
4988 (2000). 

[9] G. Giacomin and F. L. Toninelli, Phys. Rev. Lett. 96, 
070602 (2006). 



[10] F. Caravenna, G. Giacomin and M. Gubinelli, 

|math-ph/0509 065 (to be published). 
[11] A. B. Harris, J. Phys. C 7, 1671 (1974). 
[12] H. Kunz, R. Livi and A. Siito, J. Stat. Mech., P06004 

(2007). 



